Magnetic phases near the Van Hove singularity in s- and d-band Hubbard model 
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We investigate the magnetic instabilities of the nondegenerate (s-band) and a degenerate (d-band) 
Hubbard model in two dimensions using many-body effects due to the particle-particle diagrams 
and Hund's rule local correlations. The density of states and the position of Van Hove singularity 
change depending on the value of next-nearest neighbor hopping t' . The Stoner parameter is strongly 
reduced in the s-band case, and ferromagnetism survives only if electron density is small, and the 
band is almost flat at small momenta due to next-nearest neighbor hopping. In contrast, for the 
d-band case the reduction of the Stoner parameter which follows from particle-particle correlations is 
much smaller and ferromagnetism survives to a large extent. Inclusion of local spin-spin correlations 
has a limited destabilizing effect on the magnetic states. 



PACS numbers: 75.10.Lp, 71.10.-Hx, 71.28.-Hd, 75.30. Kz. 
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I. INTRODUCTION 

Although the efforts to understand the microscopic ori- 
gin of itinerant magnetism have continued for over three 
decades, there is still no consensus whether the nondegen- 
erate (s-band) Hubbard model, originally introduced as a 
model for d-electrons in transition metals,EJ can serve as a 
simple model which describes itinerant ferromagnetism. 
On the mean-field level it seems to be a good starting 
point as the Stoner criterion predicts a ferromagnetic (F) 
ground state (GS) in a broad range of parameters. How- 
ever, the inclusion joLelectron correlations destabilizes it 
in most situations.lJT3 However, the F state is stable in 
the limit of C/ ^ oo up to a criijijcal doping, being close 
to 6 = 0.29 for a square lattice. Ercl Recently, a few new 
mechanisms which stabilize ferromagnetism in the Hub- 
bard model with moderate Coulomb repulsion U have 
been proposed. They are realized either by extending 
the model by additional intersite Coulomb interactions,^ 
or in the flat-band scenario. 0^113 While the former mecha- 
nism might jwork in liquid hydrogen rather than in transi- 
tion metalsjj the question whether the spin-independent 
on-site Coulomb interaction U alone can give ferromag- 
netism remains intriguing. 

Rather extreme situations in which the GS is F are en- 
countered in low dimensional systems. First, high-spin 
GS are found in finite systems_with open-shell electronic 
states, like inr-a tetrahedron,li3 and in some other few- 
atom clusters.E^ Second, the enhanced degeneracy at the 
Fermi level realized in oupj^mension by extended hop- 
ping stabilizes tie F GS.Bil3 The mean- field analysis of 
Lin and HirschQ suggests that the F instability is en- 
hanced in two dimensions by the next-nearest neighbor 
hopping t' , and indeed the summation of the most diver- 
gent diagrams confirms the F instabiliti£_at the Van Hove 
singularity (VHS) in two dimensions .113 There are also 
indications from the enhanced stability of the Nagaoka 
state by at t/ = oo limit, that the kinetic energy chang- 
ing slowly with .electron filling favors ferromagnetism in 
a square latticeilS 



Yet, in spite of this revived interest in the magnetic 
states of the nondegenerate Hubbard model, the mag- 
netic states are realized in nature in degenerate d bands of 
3d transition metals. Therefore, we investigate here the 
F and antiferromagnetic (AF) instabilities of the d-band 
and compare them with those found in the s-band for 
the same two-dimensional (2D) lattice near the VHS. It 
is important to include the electron correlations, if mag- 
netic instabilities are considered. Using the Kadanoff- 
Baym technique of deriving conserving approximations, 
we find by including particle-particle scattering a simi- 
lar expression as-flostulated and used by Chen et al. in 
the s-band case .Ell Note that this approach is diffewjBt 
from fluctuation exchange approximation (FLEX) 
and avoids the self-consistency in the conserving approx- 
imation, but nevertheless gives a magnetic structure fac- 
tor |OfjLhc same quality as the Monte-Carlo simulatiops in 
2DJ13'E2I and in infinite dimensional Hubbard model.Ell In 
contrast, the experimental data suggest that the Hund's 
rule exchange interaction J remains practically equal to 
its atomic value,Ej but the atomic correlations stabilize 
local moments for J > 0. In the Hartree-Fock (HF) 
approximation these moments exist only in the sym- 
metry brokea phases, and they are absent in nonmag- 



netic statesEj Thus, both particle-particle scattering and 
atomic correlations contribute to the reduction of the 
magnetic energy in the d-band, conventionally expressed 
by the Stoner parameter Id- This motivated us to make a 
more extended study, in which we analyze the F and AF 
instabilities and discuss for the first time two questions: 
(i) Does ferromagnetism exist in the s-band model in a 
finite density range near the VHS? If it does, it might 
be possible to stabilize it also in three dimensions, pro- 
vided a high density of states would exist at the Fermi 
level, (ii) How does the picture change when we go to 
the d-band case? 

The paper starts with the presentation of s-band and 
d-band Hubbard models with next-nearest neighbor hop- 
ping in Sec. II. The d-band case is treated by a general- 
ization of the treatment of the particle-particle scattering 
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presented in Ref. 17. We also include atomic spin-spin 



correlations by a local ansatz method and evaluate the 
renormalized Stoner parameters. The magnetic phase di- 
agrams of the s- and d-band are presented and discussed 
in Sec. 111. The paper is concluded in Sec. IV, where we 
also give estimations of the parameters used for realistic 
transition metals. 



II. THE MODELS AND RENORMALIZED 
STONER PARAMETERS 

First, we consider a nondegenerate Hubbard (s-band) 
model on a square lattice with nearest (t > 0) and next- 
nearest neighbor {t' > 0) hopping,Ll 

(1) 

An increasing t' makes the kinetic energy, 

— — 2t(cos kx + cos ky) + At' cos k^ cos ky, (2) 

to increase slower from F to X{Y) point, and finally to 
become flat for R = 1, where R = 2t'/t. We do not 
consider the unrealistic cases with large R > 1. Second, 
we study a degenerate d-band model, with a simplified 
intraorbital hopping,cil 

(ij),ct,o- {{i,j)),a,cr 
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where = J2a'>^ia,<^ and Sia = {Sf^, Sf^, S^^} are 
density and spin operators for orbital a at site i, and 
Uq = U + 2 J. The d-bands have the same dispersion ek 
as in the s-band model (|^). 

Following the approach of Kadanoff and BaymjHj one 
may construct a 'conserving approximation' to the ri- 
band Hubbarcl— model which motivates the approach 
by Chen et aLtZl Therefore, we have considered a sys- 
tem coupled to an infinitesimal external field, hia = 



(4) 



and derived a selfenergy by taking the functional deriva- 
tive with respect to the full Green function G of the 
Kadanoff and Baym potential $, E = S^/SG. In a fi- 
nite field hia one finds the Dyson equation for the one 
particle Green's function of the d— band Hamiltonian (H), 



G (r) = G (r) - E,,, „ (r) + b, „ • a 5(r) (5) 

where G^j ^ is the noninteracting (i.e., U = J ~ 0) (di- 
agonal) Green function matrix, Sy,a is the (nondiago- 
nal) selfenergy matrix labelled by spin indices, and ct is a 
vector composed out of Pauli matrices. Due to the sym- 
metry of the hypercubic lattice considered here, there 
are no interorbital hopping processes, and therefore the 
one-particle Green function Gij^a{T) is diagonal in orbital 
space. 

The magnetic instabilities follow from the instabilities 
of the linear response function (the transverse suscepti- 
bihty) to a local spin-flip excitation {b~^ = b^^ ~ *^D' 



db-p{r) 



(6) 



b=0 



which is determined by the full one-particle Green func- 
tion G. Using a well-known identity we rewrite the func- 
tional derivative as. 



dG,,,a{T'y + e) 



b=0 



Gjii,a{T , Til) 



dbl^ir) 



, (7) 



b=0 



where / = / dr^^ stands for the integration and 

summation over all internal variables {iij}. Using the 
Dyson equation (^, we find a formal integral equation 
for the transverse susceptibility. 



b=0 



(8) 



where a = ct^ — ia"^ is a Pauli matrix, and the effective 
two particle interaction is defined as 



rap{tin,,i2TiJt3n^) = hm —— — -. (9) 

e^o+ dGl^^^j,{n„Ti,+€) 

One observes immediately, that an RPA-like expression 
for the transverse susceptibility would be obtained, if the 
selfenergy, Ey_Q,(r), were local in space and time, as is, 
e.g., the HF-selfenergy. 
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FIG. 1. Schematic representation of the diagrammatic ex- 
pansion of the Kadanoff-Baym potential part which origi- 
nates from electron-electron interaction with a representative 
d-orbital, $a- The first and second order diagrams are shown; 
higher order diagrams include multiple particle-particle scat- 
tering and dress the intraorbital Coulomb interaction Uo fur- 
ther. 



the selfenergy is obtained by summing up to infinite or- 
der the class of diagrams in the particle-particle channel, 
which leads to an alternating geometric series, shown in 
Fig. 1^. One then integrates over all internal variables 
and sums over orbitals, $ = J2a I ^a- 

Short-range electron-electron correlation effects, ne- 
glected within the RPA, consist of local particle-particle 
scattering processes. Here we use a local approximation, 
and thus the contribution of each higher order diagram 
to $Q, is oc Si J. 
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We will now use the Kadanoff-Baym technique which 
may be used to obtain conserving approximations.^ 
First, we rewrite the expression for Hd in a more 
explicit form which separates the spin-spin interaction 
into transverse and longitudinal terms, 

ia i.a<p,(T 

+ ^ [Un^a,ani/3-(T + {U - J)nia^antp^cr] 
i,a<l3.(7 

+ 4a,T^LaC,^ac,^,T- (10) 

The diagrams to lowest order in the Kadanoff-Baym func- 
tional $ (see Fig. |l|) are obtained by taking GS expecta- 
tion values of all possible contractions of the operators in 
Hd- Since we have a spin- flip term due to the infinites- 
imal field hia, also spin-flip Green functions appear. In 
choosing higher order terms, we have only kept the dia- 
grams corresponding to particle-particle scattering, sim- 
ilarly as in the FLEX method in Ref. 

For Kadanoff-Baym diagrams we only keep those which 
contribute both to the Green's function and to the sus- 
ceptibility. In this way we avoid inconsistencies when we 
calculate phase boundaries from total energy, using the 
Green's function, and from finding singularities in the 
susceptibility. This means that we only keep diagrams of 
second order in (c|^ ^c^p ^, ) , since diagrams of higher and 
lower order give vanishing contribution when we take the 
two functional derivatives and the limit of zero external 
field, hia 0. The leading diagrams of a consistent the- 
ory are shown in Fig. ^ The generating functional for 



FIG. 2. Schematic representation of the diagrammatic ex- 
pansion of the effective Kadanoff-Baym potential, with 
the effective intraorbital Coulomb interaction Uo, defined by 
Eq. ph. 



Compared to other dynamical electronic processes, 
particle-particle scattering can by considered as instan- 
taneous process. Focusing on the local particle-particle 

kernel, Xi[,a - ^-j) = GJL (n ~ ^i) (^^ " ^j), shown 
in the second-order diagrams of Fig. m, it is reasonable 



to write xT.Jn~Tj) 



where 



(11) 



with (3 = l/ksT. Within this approximation one may re- 
place the interaction between two electrons in \a j) and 
a I) states, shown by higher order diagrams in Fig. |l|, by 
the effective interaction Uq , and we can map $ onto an ef- 
fective Kadanoff-Baym potential, $ = J2aJ^a, shown 
by its diagrammatic representation in Fig. 0. The inter- 
action vertex in the effective Kadanoff-Baym potential 
obeys the following equation. 



Uo = Uo-Uox^J'iO)Uo , 



(12) 



where X^{^) = Xiiai^) ^'-'^ ^ translationally invariant 
system, which is easily solved to give an effective intraor- 
bital Coulomb interaction in d-band, 



1 + C/oxS^(0) 



(13) 



It is now straightforward to calculate the selfenergy in 
the one-particle Green's function. Performing explicitly 
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the functional derivatives, E = S^/SG, we find the local 
selfenergy matrix, 

S.,,o(r - r') = s U j Sir - r') 5,, , (14) 

with elements 

^lt = -[Uo{SrJ-jY,{S-^)]. (15) 

The adopted approximation to the selfenergy has the 
same functional form as the usual Hartree-Fock (HF) 
approximation, and therefore we call it generalized HF 
(GHF) approximation. Self-consistency for the presented 
<i>-derivable GHF approximation is required only for the 
GHF occupation numbers {riia^a)- 

Next, an RPA-like expression for the transverse sus- 
ceptibility is obtained, using the Green function with the 
selfenergy (|l|) in Eqs. (|) and (|), 

which we will call generalized (GRPA), since it has the 
same functional form as in RPA. In the s-band case the 
above formula contains Is — U (instead of Id)^ where U is 
defined by the same renormalization due to the particle- 
particle vertex as in Eq. ([l^), but with Uq replaced by U , 
and eiyes an excellent agreement with the Monte-Carlo 
dataO There seems no reason that this approximation 
should not work well also in the d-band case. 

As an important difference to RPA, the HF value of the 
Stoner parameter for d-band model, /jjf — [U -\- 6J)/5, 
is now replaced by the renormalized Stoner parameter, 

/d = i(C/o + 4J). (17) 

We note that the Stoner parameter in the s-band is just 
equivalent to the renormalized value of C/, Ig = U . 

The reduction of the Stoner parameter due to the 
screening of the intraorbital Coulomb interaction may 
be substantial in a d-band, but not quite as big as in the 
s-band case, as we show below. In contrast, the screen- 
ing of the exchange interaction J is provided by simi- 
lar expressions which involve the interorbital transitions 
on the same site, G^^^^, and is thus of second order in 

(•-L o-'-i/S o-)- If interorbital hopping vanishes (as it 
does for hypercubic lattices), J is unscreened; otherwise 
this screening is expected to be small. This is confirmed 
by the values of J deduced from the experimental data p. 
transition metals which are close to the atomic values.EHI 
We note that the renormalization of is therefore sub- 
stantially weaker than in the s-band Hubbard model dj) 
(Fig. H), where Is — U, and U is obtained from Eq. (|l2) 



with J = 0. In the 2D Hubbard model one finds that Ig is 
only weakly dependent on the band filling n, with a min- 
imum at the filling which corresponds to the VHS, but 
is finite as long as i? < 1. With increasing R, the VHS 
moves towards the lower band edge, the reduction of Ig 
gets stronger, and one finds that /s ^ for i? ^ 1. Such 
a strong renormalization of Is follows from the singular 
behavior of the particle-particle vertex. 
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FIG. 3. Stoner parameters Is/Ihf for the s-band (left), 
and Id/ Ihf for the d-band (right) as a function of filling for 
different 7i-values. The three full-drawn curves in each panel 
refer to {U + 6J)/W = 0.5, 1.0 and 2.0 from top to bot- 
tom. The dashed curves in the d-band case are results with 
spin-spin correlations (from local ansatz) included. 

The Stoner parameter is further decreased in the. 
d-band case by the atomic interorbital correlations.E^ 
This results from the formation of local moments which 
can be built into the GS by modifying aj^ngle-Slater- 
determinant HF wave function I^'hf) into,EZl 

\'^'a{{n^aA)) = cxp(-^?7,„O„0|*HF({?^»a,a})), (18) 

m 

where r]m are variational parameters, and Om are local 
operators which in the case of a d-band have the form, 

^(n) 

i,aj3 — '^ia'^ip ^ 

0\%p = S,, • S,p, (19) 

and describe local density (O-'^^^) and spin (O-^^^g) cor- 
relations, respectively. The contribution of spin correla- 
tions to the magnetic energy follows from a comparison of 
the energy obtained with the wave function \^o{{nia^a})) 
given by Ecl (|i^) with that found with density correla- 
tions only.E3 We note that the HF and correlated wave 
functions in Eq. (|l8|) are obtained for the same elec- 
tronic distribution. The Stoner parameter Id is obtained 
as a derivative of the interaction energy with respect to 
magnetization squared, and is additionally reduced by 
up to 12% of h, if {U + &J)/W = 1 (Fig. |). There- 
fore, one finds at {U + &J)/W = 1 that Id/IwF ^ 0.65 
which agrees well with the values deduced within a re- 
alistic model for 3d transition metals £3 Interestingly, at 
R — Q the largest corrections due to spin-spin correla- 
tions are found at |". — 5| ~ 2.5, which indicates that the 
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difference between local moments in the paramagnetic 
(P) and weakly F states is there larger than at n = 5.E£l 



III. MAGNETIC PHASE DIAGRAMS 

The instability of the system towards either F or AF 
order (^|) is driven by the value of the Stoner parameter 
Id (Is)- We illustrate this by considering the phases with 
a uniform and with a two-sublattice magnetic structure, 



{nia,cr) = ^ ["-0 + Km + (77 + Xai^)e'^^' 



(20) 



where m, 77 and v are order parameters, and A^- = ±1 
for a =tj i- The quasiparticles are given (up to constant 
energy shifts) by 

EL,a = ^(£k+Q + ^k) ± ^ [(£k+Q - £k)' + A^] , 

(21) 

and Q = (7r,7r) is the nesting vector at i? = 0. The 
value of the gap at half- filling (no = 1) is given by the 
effective interaction, Aa- = A = Idv/2 in the d-band case 
(Act = a = /sI//2 in the s-band case). 

The alternating magnetic order results in a two- 
sublattice magnetic structure and opens a gap Ao- for a- 
spin electrons. Assuming the filling by rio electrons (per 
one band d-subband or s-band), we analyze only the fol- 
lowing commensurate magnetic phases: (i) ferromagnetic 
(F), m = min{no, 2 — uq], v = rj = Q; (ii) partial ferro- 
magnetic (PF), m ^ 0, m < minjno, 2 — no}, v = r] — 0; 
(iii) antiferromagnetic (AF), v ^ 0, m ^ rj = 0; (iv) spe- 
cial ferrimagnetic (SFI), m = |1 — no\, u 0, tj ^ 0; here 
the stability follows from the Fermi level lying within the 
gap between two majority Slater subbands. 

In the region of their stability, the energies of magnetic 
phases are determined using the total energy expressions 
within the GHF, 



E{{nia.a}) — Eqhf + Ec. 



(22) 



where i?GHF is determined as in HF approximation from 
the quasiparticle energies (|2^). For simplicity, we give 
only the formula for less than half-filling (no < 1), 



Eg 



HF 



1 

N 



e;^ 



J2 ^ka,T 
a,keK(T) 



1 

N 



Q,keA'(i) 



(23) 



where K(^) is a set of the occupied quasiparticle states 
^ka a the lower Slater subband for cr-spin. The 

interaction energy {Hi„t) is subtracted to avoid the dou- 
ble counting, with the form of Hint determined by the 
used Hamiltonian, either or H^. 

The correlation energy depends on the magpetic order, 
and is calculated using the local ansatz [T^,ESlHZl 



In the s-band case we adopted the value of Ecorr — 0, 
which is exact for the F states, and avoids double count- 
ing of the correlation contributions in the P and AF 
states. The treatment of atomic correlations in the d- 
band beyond the GRPA is only approximate, but suffices 
to get qualitative results for the magnetic phase diagrams 
reported below. These correlations vanish in the F phase, 
but are finite in the P state, and therefore the value of the 
Stoner parameter at the F instability is found in an ap- 
proximation. On the contrary, they have no influence on 
the lines of the AF instabilities, where the order param- 
eter V increases continuosly from zero. Here we used the 
local approximation to evaluate the respective averages, 
when the exponentials in the wave functions 5'o({'T-iQ,cr}) 
( [Tsl ) are expanded in Eq. (p4|). More details may be 
found in Rcfs. |^. The magnetic phase diagrams are next 
found using the instabilities of the nonmagnetic states, 
and by comparing the total energies (^2|) of the magnetic 
phases in the region of their stability. 
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FIG. 4. Free susceptibility x''(Qi'^) as obtained at the 
nesting vector Q = (tt, tt) for the 2D densities of states with 
R = Q, 0.7, and 1.0 {W — St). The corresponding densities of 
states N{ui) are shown in the inset. 

As an illustrative example, we limit ourselves here to 
the 2D models (|^) and (||) with nearest {t) and next- 
nearest (t') neighbor hopping. By changing the value of 
R = 2t'/t, the bands (p^ become flat and the noninter- 
acting susceptibihties x^(k, w) change. The singularity 
in the AF susceptibility x"(Q,^^) [Q = (Tr, tt)] moves to 
lower energies, as shown in Fig. |, and gradually dis- 
appears, while the FM susceptibility, x'^{0,(^) — N{u!), 
develops a peak at low energies (see the inset of Fig. ^) . 
Therefore, the RPA instabilities of the s-band model to- 
wards FM states occur at all electron densities (Fig. ||), 
and are enhanced at low n for increasing t'.u One finds 
that these instabilities occur towards the saturated FM 
states (with a maximum value of m = |1 — no|) in most 
cases. As expected, the AF (SFI) order is more stable 
near half-filling. 
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If correlation effects (renormalization of U within the 
GRPA) are included, the phase diagrams change dras- 
tically (Fig. H). The FM instability disappears almost 
entirely in the GRPA if i? < 1, except just in a narrow 
region around the VHS for larger values of R. At i? = 
it is suppressed at any n, in agreement with Rudin and 
Mattis£j As also no instability of the P ^tatc towards the 
FM order is found in infinite dimension£3l it is likedtt that 
it does not exist in hypercubic lattices at i? = O.Eil The 
situation changes when particle-hole symmetry is broken 
for t' > 0, and the FM order is stabiiiaed by infinitesi- 
mal U at the VHS point, if i? > 0.55Jl3 Remarkably, at 
R = 1 ferromagnetism is stable at low density, n < 0.45, 
as a consequence of flat band behavior, £k = — 2f, if 
k= (fc^,0) ork= (0,fcy). 




FIG. 5. RPA (left) and GRPA (right) phase diagrams for 
the s-band case and for R = 0, 0.4, 0.7, and 1.0 from top to 
bottom. FuU lines give the instabilities with increasing U of 
paramagnetic (P) to ferromagnetic (F) or antiferromagnetic 
(AF) states. The special ferrimagnetic states (SFI) are sepa- 
rated by dot-dashed lines. Ferromagnetic states are saturated 
(F) except for small regions in the three lowest panels to the 
left (PF). 

The AF order is found to be more robust and exists 
around half- filling at any value of i? < 1, with elec- 
tron correlations included. In agreement with expec- 
tations, the region of AF order at intermediate values 
of U/W > 1.0 expands when the kinetic energy in the 
ordered state increases with the increasing intrasublat- 
tice hopping both in RPA and in GRPA, as shown 
in Fig. |5[ On the contrary, at i? = 1 the region 
of antiferromagnetism is much reduced compared with 
R = 0. It is expected that incommensurate magnetic 
order with AF correlations is stable in between the AF 
and P phases, as shccsin recently within this method in 
infinite dimension OS 

In contrast to ferromagnetism, the tendency towards 
AF order at and near half-filling is weakened by the in- 
creasing next-nearest neighbor hopping t' , as the shape 



of the Fermi surface changes and the perfect nesting con- 
dition is not satisfied. At i? = 0, the transition to the 
AF state occurs at infinitesimal U > 0, and the magnetic 
moment v increases gradually from zero with increasing 
U. This behavior is characteristic of perfect nested band 
structures, and is replaced by a jump to a finite magne- 
tization v which occurs at finite [/ > 0, if i? > 0. With 
increasing R the lower quasiparticle subband E~ (k) ( |2l| ) 
becomes more flat and it is gradually more difficult to 
stabilize an AF state. Therefore, the critical value of the 
Coulomb interaction Uc/W, increases with increasing R, 
as shown in Fig. ^ 
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FIG. 6. Instability of a paramagnetic metallic (PM) to an- 
tiferromagnetic insulating (AFI) ground state at half-filling 
(n = 1) in the Hubbard model (1): (a) critical value of the 
Coulomb interaction Uc as function of i? = 2f'/ii (b) elec- 
tronic structure in an almost insulating (weakly metallic) AF 
state obtained at (7 = 4.04VF and R ^ 1. 

It is interesting to ask whether the insulating AF state 
at half filling and i? = is replaced by itinerant anti- 
ferromagnetism with increasing values of t' . The gap in 
the AF state at the transition point is determined by the 
energy difference between the top of the lower Slater sub- 
band E" at k = (7i'/2, 7r/2), and the minimum of the up- 
per Slater subband at k = (tt, 0). Using the effective 
Coulomb interaction U, we find that the AF state is al- 
ways insulating at the transition point, except at i? = 1, 
where the gap between the Slater subbands closes and 
the AF state is weakly metallic, as shown in Fig. 0. The 
onset of metallic behavior occurs in GRPA at 2t = A, 
where = A = Isv/2_m Eq. (|l[). If i? = 1, one finds 
at the transition point Uc = 0.409VF, and v = 0.606, and 
the splitting A = 0.124VF, just somewhat smaller than 
2t' — 0.125VF. This results in a very small difference of 
the energies, AE = E+{X) - E-{S) = -2.2 x lO^^tf, 
and supports the recent conclusion of Duffy and MoreoEHl 
that itinerant magnetism is difhcult to realize in the 
U — t — t' Hubbard model, without the second neigh- 
bor hopping t along the x and y directions. Morever, 
it seems that the critical values of U obtained for the 
onset of AF LRO might be overestimated at finite R, 
as we find the AF instability at Uc — 4.5t for R = 0.4, 
while the quantum Monte-Carlo calculations suggest a 
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value Uc — 2.5t. IS However, they agree with the metal- 
insulator transition estimated to be somewhere in the 
range of U < Uc < 6t, if i? = 0.4. 

The RPA phase diagram of the d-band model is simi- 
lar to that of the s-band (Fig. |^) , but the SFI phase is 
destabilized by the interorbital exchange interaction J, 
and thus the AF state is found instead in a broader range 
of parameters. As a consequence of the weaker screen- 
ing by particle-particle diagraais the conditions for 
FM LRO are less restrictive)^ and the FM instabilities 
occur in the GRPA at any value of R, but for larger in- 
teractions (Fig. 1^. As in the s-band case, FM order is 
favored at low electronic filling for 0.4 < R < 1.0 due 
to weaker dispersion. As already seen in the screening 
of the Stoner parameter Id, the spin-spin correlations re- 
strict the region of FM states, in particular at i? = 0, 
and for n > 7.5, if R > 0.4. Furthermore, one finds that 
weak ferromagnetism is somewhat more pronounced than 
in the s-band and survives for the screened interactions, 
but occurs only in a relatively narrow range of n. This 
shows that partly polarized FM states are more likely to 
result from either the local maxima or the splittings be- 
tween Cg and t2g orbitals in realistic band structures of 
3d transition metals. 
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FIG. 7. The same as in Fig. ^, but for the d-band model 
(2) with J/U — 0.25. The magnetic instabilities are shown 
bu full lines, while dotted lines indicate the same instabilities 
if the spin-spin correlations are neglected. 



IV. SUMMARY AND CONCLUSIONS 

Summarizing, we conclude that the FM states are sta- 
ble not only at, but also in the neighborhood of the VHS 
in the s-band model, if the band filling n is small. This 
suggests that ferromagnetism might be realized in the 
Hubbard model also in three dimensions, but only if the 
particle-hole symmetry is broken and the kinetic energy 
satisfies rather extreme conditions. Therefore, the s-band 
model cannot serve as a generic description of itinerant 



magnetism in transition metals. Instead the orbital de- 
generacy of the d-states is crucial to explain magnetic 
states, and the interorbital exchange coupling J plays an 
important role. 

Although no more than qualitative statements can be 
made for realistic transition metals, it is interesting to 
compare the obtained phase diagrams with the known 
interaction parameters of 3d systems. The values of U 
and J are approximately known for 3d transition met- 
als and may be obtained from the multiplet splittings, 
as discussed in detail by van der Marel and Sawatzky.E3 
This analysis leads to U = 2.55, 2.76, and 2.97 eV, for 
Fe, Co, and Ni, respectively, while tbe values of J are 
given by the same ratio, J/U = 0.27.c3 These values of 
U are close to those considered by one ot-us in an ear- 
lier study of ferromagnetism in 3d metals,E3 and we treat 
them here as representative ones. Taking the bandwidths 
of Fe, Co-and Ni as determined by Andersen, Jepsen, and 
GlotzehEl W = 5.43, 4.84, and 4.35 eV, respectively, this 
results in the HF values of the magnetic interaction for 
these elements of {U + 6J)/W = 1.23, 1.49, and 1.79, 
respectively. We assumed the electron filling of n = 7.3, 
8.3, and 9.4 in the d-band and i? = in our model ^ 
to simulate qualitatively the situation in Fe, Co, and Ni, 
and found a reduction of the respective Stoner parameter 
due to particle-particle renormalization and the spin-spin 
correlations of the order of Id/lHF — 0.61, 0.58, and 0.54, 
respectively. Interestingly, these values are almost con- 
stant as the increasing values of {U + 6J)/W are counter- 
balanced by the smaller reduction of Id (both by particle- 
particle and even more by spin-spin correlations) when 
the electronic filling increases towards the filled band. 
Taking JJII = 0.25 in a realistic model with canonical 
d bands,tZl one finds instead that Id/IuF — 0.63, 0.57, 
and 0.52. The agreement is very good indeed granted 
the simplicity of the model. 

Furthermore, using the values |Of Id derived above and 
realistic values of d-bandwidthsj2a one finds the Stoner 
parameters: Id ^ 0.72, 0.82, and 0.92 eV for Fe, Co, 
and Ni, respectively. These values are significantly lower 
than those deduced by Gunnarsson from the local den- 
sity approximation (LliA) calculations, being 0.92, 0.99, 
and 1.01, respectively.L^I In spite of the qualitative na- 
ture of this comparison, this allows us to conclude that 
important corrections of the Stoner parameter exist due 
to nonlocal electron correlation effects, as in particular 
due to the spin-spin interorbital correlations, which can- 
not be dealt with in the standard band structure calcula- 
tions based on the LDA. Thus, in spite of some attempts 
which exist in the literaturej2Zl there is no way to deduce 
reliable values of the Stoner parameter from band struc- 
ture calculations performed within the LDA. 

An additional factor which might stabilize ferromag- 
netism in the d-band model is the fiattening of the bands 
with increasing values of R. Then magnetic states are 
possible even for rather small interactions U. Of course, 
this tendency is overemphasized in the present 2D model 
by the Van Hove singularity, but it is expected that 
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strong next-nearest neighbor hopping t' may lead to FM 
instabihties also in three-dimensional or quasi 2D sys- 
tems at low filling. It would be interesting to verify this 
prediction, if such materials could be synthesized. 

Altogether, we have shown that a correct quantita- 
tive description of ferromagnetism in transition metals 
is only possible within a realistic d-band model and 
when the particle-particle screening and spin-spin corre- 
lations are included. In order to obtain more quantitative 
results, however, realistic densities of states have to be 
used. A somewhat different situation, however, is found 
for the AF states, where the band structure effects (like 
perfect nesting) dominate, at least at weak and inter- 
mediate values of U , and the s-band model might then 
suffice to explain quajlitatively the experimental data for 
AF Cr and its alloys EI 
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